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Abstract 

We show that flavor effects in leptogenesis reduce the region of the see-saw 
parameter space where the final predictions do not depend on the initial conditions, 
the strong wash-out regime. In this case the lowest bounds holding on the lightest 
right-handed (RH) neutrino mass and on the reheating temperature for hierarchical 
heavy neutrinos do not get relaxed compared to the usual ones in the one-flavor 
approximation, Mi(T reh ) > 3(1.5) x 10 9 GeV. Flavor effects can however relax 
down to these minimal values the lower bounds holding for fixed large values of the 
decay parameter K%. We discuss a relevant definite example showing that, when the 
known information on the neutrino mixing matrix is employed, the lower bounds for 
K\ 10 are relaxed by a factor 2-3 for fully hierarchical light neutrinos, without 
any dependence on #13 and on possible phases. On the other hand, going beyond 
the limit of hierarchical light neutrinos and taking into account Majorana phases, 
the lower bounds can be relaxed by one order of magnitude. Therefore, Majorana 
phases can play an important role in leptogenesis when flavor effects are included. 



1 Introduction 



The see-saw mechanism [TJ is a simple and elegant way to understand neutrino masses 
and their lightness compared to all other fermions. At the same time, with leptogenesis 
[2], the see-saw offers an attractive model to explain the matter-antimatter asymmetry of 
the Universe. The asymmetry is produced in the decays of the very heavy RH neutrinos, 
necessary for the see-saw mechanism to work, generating a B — L asymmetry in the form 
of a lepton asymmetry partly converted into a baryon asymmetry by sphaleron processes 
if temperatures are larger than about 100 GeV [3]. 

Calculations have been typically performed within a one-flavor approximation, assum- 
ing that the final total B — L asymmetry is not sensitive to the dynamics of the individual 
B/3 — L a asymmetries. In this case it has been noticed [HE] that current neutrino mixing 
data favor leptogenesis to lie in a region of the parameter space where thermal equilibrium 
sets up prior to the freeze-out of the final asymmetry and therefore, if the temperature 
of the early Universe is large enough, the final asymmetry does not depend on the initial 
conditions. 

For a hierarchical heavy neutrino spectrum, successful leptogenesis requires the mass 
of the lightest [6] or of the next-to-lightest [7] RH neutrino to be large enough, if the RH 
neutrino production occurs through scatterings in the thermal plasma. In this case there 
is also an associated lower bound on the initial temperature of the radiation-dominated 
regime in the early Universe, to be identified with the reheating temperature, T re h, within 
inflation. Being independent of the initial conditions, these bounds are particularly rele- 
vant, since they appear as an intrinsic feature of leptogenesis. 

Recently it has been pointed out that flavor effects can strongly modify leptogenesis 
predictions 0, EE Ell El 021 E31 [Ej. In this work we investigate how the dependence 
on the initial conditions and the lower bounds on the RH neutrino masses and on the 
reheating temperature are modified when flavor effects are included. We show that the 
usually quoted values do not get relaxed, contrarily to what has been stated in [13] . The 
calculations are presented within the hierarchical limit for the spectrum of the heavy 
neutrino masses, and more specifically in the lightest RH neutrino dominated scenario 
(iViDS), where the final asymmetry is dominantly produced by the lightest RH neutrino 
decays. In the second Section we set up the notation and the kinetic equations. In the 
third Section we analyze how flavor effects modify the dependence on the initial conditions. 
In the fourth Section we perform explicit calculations for a specific but significant example. 
In the Appendix we discuss the role of AL = 1 scatterings, showing how they introduce 
just a correction to the results. 
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2 General set-up 



Adding to the Standard Model Lagrangian three RH neutrinos with a Yukawa coupling 
matrix h and a Majorana mass matrix M, a neutrino Dirac mass matrix mp = hv is 
generated, after electro-weak symmetry breaking, by the vev v of the Higgs boson. For 
M ^> m D , the neutrino mass spectrum splits into 3 heavy Majorana states Ni, N 2 and 
A/3 with masses Mi < M 2 < M3, almost coinciding with the eigenvalues of M, and 3 light 
Majorana states with masses mi < m 2 < m^, corresponding to the eigenvalues of the 
neutrino mass matrix given by the see-saw formula, 

m v = -m D —m T D . (1) 

Neutrino mixing experiments measure two mass-squared differences. In normal (inverted) 
neutrino schemes, one has 

m\-m\ = Am 2 atm {Am 2 sol ) , (2) 

m\-m\ = Am 2 sol (Am^ tm - Am 2 sol ) . (3) 

For mi ^> m atm = a/ Am 2 tm + Am 2 ol ~ 0.05 eV, one has a quasi-degenerate spectrum 
with mi ~ m 2 — m 3 , whereas for mi <C m so \ = a/ Am s 2 ol ~ 0.009 eV one has a fully 
hierarchical (normal or inverted) spectrum. 

A lepton asymmetry can be generated from the decays of the heavy neutrinos into 
leptons and Higgs bosons and partly converted into a baryon asymmetry by the sphaleron 
(B — L conserving) processes at temperatures higher than about 100 GeV. 

For 10 9 GeV < Mi < 10 12 GeV, tauon charged lepton Yukawa interactions are in 
equilibrium [HI [13l [12] and typically faster than inverse processes [15], breaking the 
coherent evolution of the lepton doublets quantum states, the |£j)'s, produced in RH neu- 
trino decays. The quantum state is projected on a two-flavor basis, where the eigenstates 
are given by the tauon flavor and by a linear combination of muon and electron flavors. 
For Mi < 10 9 GeV both muon and tauon charged lepton Yukawa interactions are strong 
enough that the quantum state is projected on a three-flavor basis. 

Flavor effects play a double role in leptogenesis. A first effect is that the wash-out 
is reduced. This happens because, in the inverse decays, the Higgs do not interact with 
|/j)'s but with the flavor eigenstates |/ a )'s, with a reduced inverse decay rate. This effect 
can be accounted for introducing the flavor projectors [H [12] 

Pia = |^|U| 2 = ^, (4) 
1 See also the note added in the end of the paper. 
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P ia = K/:,|L)| 2 = ^, (i = l,2,3;a = e,/i,r) (5) 

where r iQ , is the partial decay rate of the process N { — > l a + W and Tj = J2 a is the 
total decay rate, such that ^ Q P iQ = 1. Analogously T ia is the partial decay rate of the 
process N { — > l a + if with Tj = ^ Q r ia , such that J2 a ^ = 1- 

A rigorous description of the asymmetry evolution has then to be performed in terms 
of the individual flavor asymmetries A a = B/3 — L a rather than in terms of the total 
asymmetry Nb-l = Yla -^A a) as usually done in the one-flavor approximation. The con- 
tribution of each decay to the N&Js is determined by the individual flavor CP asymmetry, 
defined as 

F. — P. 

B ia = ~= ■ (6) 

A second role played by flavor effects arises because the state is not the CP conjugated 
state of \k) [HI [12] and this yields an additional source of CP violation. This can be 
described in terms of the projector differences [SJ E2] 

AP = P — P (7) 

L - 1± ia — 1 ia 1 ia i \ ' J 

obeying £ a AP ia = 0. Indeed writing P ia = P° Q + AP ia /2 and P ia = P° Q - AP ia /2, 
where the Pf a = (P ia + P ia )/2 are the tree level contributions to the projectors, one can 
see that 

A P ia 

e ia = eiP ia -\ - — , (8) 

where the total CP asymmetries £j = ^ a Si a . The three flavored CP asymmetries can 
be calculated using [16] 



(9) 



where Xi = (Mi /Mi) 2 and 



£(x) = ^x 



'1 + x) In 



1 + x \ 2 — x 



x J 1 — X 



(10) 



In the last years there has been an intense study of the relevant processes in leptogenesis 
[T7J [8j [18j [TH [20]. For our purposes it will be enough to describe the evolution of the 
asymmetries just in terms of decays and inverse decays (after subtraction of the real 
intermediate state contribution from AL = 2 processes) neglecting AL = 1 scatterings 
and thermal effects. In the Appendix we will discuss and justify this approximation. 
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We will also neglect off-shell AL = 2 and AL = processes, relevant only at higher 
temperatures, and spectator processes that produce corrections not larger than 30%, 
within the precision needed for our purposes [20]. The kinetic equations, using z = M\/T 
as the independent variable, are then given by [HI CGI EH [21] 

= -DiiNn-Nfi) (11) 

i i 

where the N&Js and the iVj's are the abundances per number of iVx's in ultra-relativistic 
thermal equilibrium. Notice that a = e, fi, r in the three-flavor case, while a = r, e + \i 
in the two-flavor case, where 'e + /i' means that the electron and muon CP asymmetries 
Ei a and projectors Pf a , are summed. Therefore, in the two-flavor case, there are only 
two kinetic equations for the N& a instead of three. The equilibrium abundances are 
given by iV^ = zf /C 2 (^)/2, where we indicate with /Q(zj) the modified Bessel functions. 
Introducing the decay parameters K { = T D i {T = 0)/H(T = Mi), the ratios of the total 
decay widths to the expansion rate at T = Mj, the decay terms Di can be written like [5] 

D t = ^=K l x l z (-) , (12) 



Hz \7, 

where T Dji = r, + Tj = T Dji (T = 0) (l/"fi) are the total decay rates and the (l/7j)'s are 
the thermally averaged dilation factors and are given by the ratios K\{zi)jKi{zi). Finally, 
the inverse decays wash-out terms are given by 

Wl D (z) = ±K l ^lC 1 (z l )zl (13) 

The evolution of the N^s can be worked out in an integral form, 

N Aa (z) = P - ^ dz ' W ^ + Y, ^ *«(*) , (14) 

i 

with the 6, in the two-flavor case, or 9, in the three-flavor case, efficiency factors given by 
(*; K it Pi) = - f dz' ^ e-a ^ . (15) 



The total final B — L asymmetry is then given by N B _ L = . Finally from 

this, assuming a standard thermal history and accounting for the sphaleron converting 
coefficient a sp h ~ 1/3, the final baryon-to-photon number ratio can be calculated as 

Vb = a sph * 0.96 x 10- 2 N B _ L , (16) 
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to be compared with the measured value [22J 

r^MB = ( 6 .3 ± 0.3) x 1(T 10 . (17) 

Notice that the efficiency factors depend only on the Pf a but not on the differences AP ia . 
The final asymmetry depends in general on all the unmeasured 14 see-saw parameters. A 
useful parametrization is provided by the orthogonal see-saw matrix Q [23J, in terms of 



which one has mjj = U \/D m Q V 'Dm- In this way the decay parameters can be expressed 
as 

i 

where m+ ~ 10 -3 eV is the equilibrium neutrino mass, and the tree level projectors as 



It will also prove useful to introduce the quantities K ia = P® a Ki. The orthogonal matrix 
can in turn be decomposed as 

ft = R12M R 13 (u 31 ) R 23 (u 22 ) , (20) 

such to be parameterized in terms of three complex numbers, u 2 i, ou^i and u 22 , deter- 
mining the rotations in the planes 12, 13 and 23 respectively. In this way one has 
Vb = VB( m i,U, Mi,u>ij). It is interesting that including flavor effects there is a poten- 
tial dependence of the final asymmetry also on the unknown parameters contained in the 
PMNS mixing matrix U [9lfT2]. 

n the following we will assume a hierarchical heavy neutrino spectrum with M 3 , M 2 
In this case the general expression Eq. (J9j) for the CP asymmetries E\ a , reduces to 



M 



8 7T(/lU) n 



Er 



3 1 
h* al h ui ( ——{tih) Xj + —{h)h) jX 



(21) 



Just in passing, we notice that it is straightforward to generalize, including flavor effects, a result 
obtained in [5T] for the efficiency factors in the degenerate limit, obtained for (M3 — Mi) /Mi < 0.1, 
within the one-flavor approximation and for Ki 3> 1, for all i. Indeed, in this case one can approximate 
dNi/dz 1 ~ dN^/dz' in the Eq. (fT"5)) obtaining that n\ a ~ k(K ), where K a = J2i Kia- The function 
k{x) is defined in the Eq. (|27|) and it approximates n\ a in the hierarchical limit when x — Ki a . In the 
degenerate limit one has then just to replace Ki a with the sum K a . Like in the hierarchical limit, the 
number of efficiency factors to be calculated reduces from 6 to 2 in the two-flavor case and from 9 to 3 in 
the three-flavor case, one for each flavor. If instead of a full degeneracy, one has just a partial degeneracy, 
with M 3 3> M 2 — Mi, then simply one has K f 3a <C n la ~ n l 2a ~ n(Ki a + K 2a )- 
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We will moreover assume no rotation in the plane 23, i.e. R23 = 1. Under these conditions 
both the total CP asymmetries e 2 and e 3 are suppressed ~ Mi/M 2 ^. It is interesting to 
notice that, under particular conditions, the £2a's and the e^a's are not suppressed and 
do not vanish. This can potentially lead to a scenario where the final asymmetry is 
produced by the decays of the two heavier RH neutrinos if M 2 3 < 10 12 GeV such that 
the flavored regime applies. Here we do not pursue this possibility and focus on a typical 
A^i-dominated scenario where the dominant contribution to the final asymmetry comes 
from the decays of the lightest RH neutrino and 



N b-l - 2. £ ^ < ■ (22) 

a 

It will prove important for our discussion that the CP asymmetries, the total e% and 
the flavored £i a 's, cannot be arbitrarily large. The total CP asymmetry is indeed upper 
bounded by El S] 

ex <e(M 1 )(3 max f(m 1 ,K 1 ), (23) 

where e(Mi) = 3 M x m atm /(16 ix v 2 ), /3 max = m atm /(mi + m 3 ) and < /(mi, K{) < 1 is a 
function that vanishes for K\ = m 1 /m* and tends to unity for mi/ (m* Ki) — > 0. On the 
other hand each individual flavor asymmetry e la is bounded by [1] 



|£ lQ | < e(Mi) JP^-^-maXj [\U aj \] (24) 

V m atm 

and therefore while the total CP asymmetry is suppressed when mi increases, the single 
flavor asymmetries can be enhanced. The existence of an upper bound on the quantity 
ri a = £i a /e(Mi) independent of Mi, imply, as in the case of one-flavor approximation, 
the existence of a lower bound on Mi given by 

Ml > Mr ,»^Hp!^W >!^i^ (at3 .C.L.), ,25) 
a sph 3 k\(Ki) ^(Ki) k\(Ki) £i(Ki) 

where we indicated with k\(Ki) the efficiency factor in the usual one-flavor approximation, 
corresponding to n\ a with P° Q = 1. We also defined 

6 = with 6^ ria f;^ la) - (26) 

This quantity represents the deviation introduced by flavor effects compared to the one- 
flavor approximation in the hierarchical light neutrino case. We could then use ri a < 
\/ Pi a m 3 1 m atm to maximize ^1. Notice however that first, because of the bound on the 
total asymmetry, the r la 's cannot be simultaneously equal to \fP~i a and second in ^ there 



are sign cancellations. Therefore, the bound Eq. (125]) is more restrictive than this possible 
estimation and we prefer to keep it in this form, specifying £1 in particular situations. 

As usual, the lower bound on Mi imply also an associated lower bound on the reheating 
temperature T re h. 



3 Dependence on the initial conditions 



From the Eq. (115j) . extending an analytic procedure derived within the one-flavor ap- 
proximation [5], one can obtain simple expressions for the k( q 's. In the case of an initial 
thermal abundance (N^ = 1) one has 

K\ a ~K{K Xa )= (l-e ^ ), (27) 

Ki a ZB{Ki a ) V / 

where 

^ ltt )-2 + 4Ce"^. (28) 

Notice that in the particularly relevant range 5 < K\ a < 100 this expression is well 
approximated by a power law [25J, 

0.5 



in , vl -2 (29) 

In the case of initial vanishing abundance (N^ = 0) one has to take into account two 
different contributions, a negative and a positive one, 

K\ a = d_(K u Pi) + ^(K,, Pi) . (30) 

The negative contribution arises from a first stage when < Nff , for z < 2^, and is 
given approximately by 

Kl( Kll Pl) * -A _ A . (31) 

The iVj abundance at z eq is well approximated by 

N Nl ( Zeq ) ~ = - , (32) 

interpolating between the limit K\ ^> 1, where z eq <C 1 and N^ 1 (z eq ) = 1, and the limit 
K\ <C 1, where z eq ^> 1 and NN 1 (z eq ) = N(Ki) = 37riCi/4. The positive contribution 
arises from a second stage when iVjvi > f° r £ > ^eq, and is approximately given by 

2 / K la z B (K la )N N (z ccl )\ 

K f + (^ 1; J* ) ~ - 1 - e ■ . (33) 

z B {K la )Ki a \ J 
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Notice, from the Eq. (1321) . how Njsr^Zeq) is still regulated by K\ and this because the RH 
neutrino production is not affected by flavor effects, contrarily to the wash-out, which is 
reduced and is regulated by K\ a . 

These analytic expressions make transparent the two conditions to have independence 
on the initial conditions. The first is the thermalization of the N\ abundance, such 
that, for an arbitrary initial abundance one has N^^Zeq) = 1. The second is that the 
asymmetry produced during the non-thermal stage, for z < z eq , has to be efficiently 
washed out such that <C They are both realized for large values of Ki ^> f. 
More quantitatively, it is useful to introduce a value K ± such that, for K\ > K*, the final 
asymmetry calculated for initial thermal abundance differs from that one calculated for 
vanishing initial abundance by less than some quantity 5. This can be used as a precise 
definition of the strong wash-out regime. Let us consider some particular cases, showing 
how flavor effects tend to enlarge the domain of the weak wash-out at the expense of the 
strong wash-out regime. 

3.1 Alignment 

The simplest situation is the alignment case, realized when the iVi's decay just into one- 
flavor a, such that P la = Pi a = 1 and P\p^ a = Pip^ a = 0, implying e la = e\. Notice that 
we do not have to worry about the fact that the lightest RH neutrino inverse decays might 
not be able to wash-out the asymmetry generated from the decays of the two heavier ones, 
since we are assuming negligible £2/3 and anyway. In this case the general set of kinetic 
equations (ITT]) reduces to the usual one-flavor case and all results coincide with those in 
the one-flavor approximation [12] . In particular one has N { B _ L = e\ K la . 

In the considered case of alignment, one obtains K ± ~ 3.5 for 5 = 0.1, as shown in 
Fig. 1. The value of plays a relevant role since only for K\ > one has leptogenesis 
predictions on the final baryon asymmetry resulting from a self-contained set of assump- 
tions. On the other hand for K\ < leptogenesis has to be complemented with a model 
for the initial conditions. 

In addition to that, one has to say that in the weak wash-out regime the calculation 
of the final asymmetry requires a precise description of the RH neutrino production, 
potentially sensitive to many poorly known effects. It is then interesting that current 
neutrino mixing data favor K\ to be in the range K so \ ~ 9 < K\ < 50 ~ -Katm [H EE U\, 
where one can have a mild wash-out assuring full independence of the initial condition, as 
one can see in Fig. 1, but still successful leptogenesis. In this case one can place constraints 
on the see-saw parameters not depending on specific assumptions for the initial conditions 



9 



10"' 10° 10' 10 2 10 3 10" 10° 10' 10 2 10 3 




K 1 K , 

Figure 1: Efficiency factor (left panel) and lower bounds on M\ (right panel) in the 
alignment (solid), in the semi-democratic (dashed lines) and in the democratic (short- 
dashed) cases. 

and with reduced theoretical uncertainties. 

Since one has £i = 1, the general lower bound on M\ (see Eq. ff25l) ). like all other 
quantities, becomes the usual lower bound holding in the one-flavor approximation. In 
the right panel of Fig. 1 we have plotted it both for = 1 and = 0. One can see 
that the K la dependence on the initial conditions translates into a dependence of M™ in . 
The lowest model-independent values are then obtained for K\ = K+ ~ 3.5 and are given 
by 

M 1 > 3 x 10 9 GeV and T reh > 1.5 x 10 9 GeV . (34) 

Another typically quoted lower bound on Mi is obtained for initial thermal abundance in 
the limit K 1 -> 0, given by Mi > 4 x 10 s GeV [2S]. 

3.2 Democratic and semi-democratic cases 

Let us now discuss another possibility and for definiteness let us first consider the three- 
flavor regime. We will easily extend the results to the two-flavor regime. Let us assume 
that Pi a = Pi a = 1/3 for any a and consequently APi a = 0. This case was also 
considered in [13]. From the Eq.'s (fl5|) it follows that the three n\ a , like also the three 
Eicxi are all equal and thus the Eq. (|22j) simplifies into N { B _ L = ei K\ a , as in the usual 
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one-flavor approximation. However now the wash-out is reduced by the presence of the 
projector, such that K% — > P® a Ki = Ki/3. The result is that in the case of an initial 
thermal abundance the efficiency factor, as a function of K\, is simply shifted. The same 
happens for vanishing initial abundance in the strong wash-out regime. However, in the 
weak wash-out regime, there is not a simple shift, since the RH neutrino production is 
still depending on K\. A plot of n\ a is shown in the left panel of Fig. 1 (dashed lines). 
One can see how the reduced wash-out increases the value of to ~ 10, approximately 
1/Pi a — 3 larger, thus compensating almost exactly the wash-out reduction by a factor 
~ 3 1 ' 2 (cf. (J29D). 

In this way the lowest bounds on Mi in the strong wash-out regime, at Ki = K+, are 
almost unchanged. On the other hand for a given value of Ki ^> K+, the lower bounds 
get approximately relaxed by a factor 3 |13j . The lower bounds for the democratic case 
are shown in the right panel of Fig. 1 (dashed lines). Notice that also the lower bound 
at Ki —>■ for initial thermal abundance does not change. More generally, flavor effects 
simply induce a shift of the dependence of the lower bound on Ki. 

The semi-democratic case is intermediate between the democratic and the alignment 
cases. It is obtained when one projector vanishes, for example Pip = 0, and the other two 
are one half and in this case K* ~ 7. The corresponding plot of the efficiency factor and 
of the lower bound on Mi are also shown in Fig. 1 (short-dashed lines). These results 
for the semi-democratic case also apply to the two-flavor regime, when the two projectors 
are equal, namely P° = -P e ° +A1 = 1/2. 

3.3 One-flavor dominance 

There is another potentially interesting situation that motivates an extension of the pre- 
vious results to arbitrarily small values of P° a . This occurs when the final asymmetry is 
dominated by one flavor a and the Eq. (1221) can be further simplified into 

N' B _ L ~ e lQ n[ a , (35) 

analogously to the alignment case but with P® a <C 1. Notice that this cannot happen due 
to a dominance of one of the CP asymmetries, for example with ei a being close to its 
maximum value, Eq. ( J2^|) . much larger than the other two strongly suppressed, simply 
because one has Y2 a A-Plo = 0. One has then to imagine a situation where the CP 
asymmetry e lQ is comparable to the sum of the other two but K lf s^ a ^> K la > 1 such 
that K la ^> k\q. The dominance is then a result of the much weaker wash-out. 

The analysis of the dependence on the initial conditions can then be again performed 
as in the previous cases calculating the value of for any value of Pf a . The result is 
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Figure 2: Value of K+, defining the strong wash-out regime for 5 = 10% (solid line) and 
5 = 50% (dashed line). 

shown in Fig. 2. The alignment case corresponds to P® a = 1, the semi-democratic case to 
P\ a = 1/2 and the democratic case to P° a = 1/3. Notice that the result is very close to 
the simple estimation K*(P® a ) = K+(l) / P® a , that would follow if n\ a were just depending 
on K la . In Fig. 3 we have plotted the values of the lower bounds on Mi and T re h for 
light hierarchical neutrinos, implying m 3 = m atm in the Eq. (jUj) . These can be obtained 
plugging £i = \/Ffa K\ a (K\ a ) / k\(Ki) < 1 in the Eq. ( |25i) . They correspond to the lowest 
values in the strong wash-out regime, when K\ > K±. 

There are two possible ways to look at the results. As a function of P° a they get more 
stringent when P® a decreases so that the minimum is obtained in the alignment case, 
corresponding to the unflavored case. This is clearly visible in the left panel of Fig. 3. 
Therefore, one can conclude that flavor effects cannot help to alleviate the conflict of the 
leptogenesis lower bound on the reheating temperature with the gravitino problem upper 
bound. On the other hand, flavor effects can relax the lower bound for fixed values of 
K\ 1. Indeed, for each value K\ 3> 1, one can choose P? a = K*(P? a = l)/K 1 such that 
Ki a = K*(Pia)i ^us obtaining the highest possible relaxation in the strong wash-out 
regime. The result is shown in the right panel of Fig. 3. It is important to say that this 
relaxation is potential. A direct inspection is indeed necessary to determine whether it is 
really possible to achieve at the same time not only small values of P® a , but also a £\ a 
that is not suppressed compared to E\p± a . 

Notice that together with the one-flavor dominance case, one can also envisage, in 
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Figure 3: Lower bounds on Mi and T re h calculated choosing P® a = K+(P® a = 1)/Ki such 
that K\ = K^Pia) (thick lines) and compared with the usual bounds for P° a = 1 (thin 
lines). In the left panel they are plotted as a function of P° a , while in the right panel as 
a function of K\. 

the three-flavor regime, a two- flavor dominance case, where two projectors are equally 
small, with the third necessarily close to one, while all the 3 flavored CP asymmetries are 
comparable. 

In the next Section, we will consider a specific example that will illustrate what we 
have discussed on general grounds. At the same time it will help to understand which 
are realistic values for the projectors and their differences, given a specific set of see-saw 
parameters and using the information on neutrino mixing matrix we have from low-energy 
experiments. 

4 A specific example 

The previous results have been obtained assuming no restrictions on the projectors. More- 
over, in the one-flavor dominance case, where there can be a relevant relaxation of the 
usual lower bounds holding in the one-flavor approximation, we have assumed that the 
upper bound on e la , see Eq. (J24l) . is saturated independently of the value of the projector. 

This assumption does not take into account that the values of the projectors depend 
on the different see-saw parameters, in particular on the neutrino mixing parameters, 
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and that severe restrictions could apply. Let us show a definite example considering a 
particular case for the orthogonal matrix, 



n = r 



13 



/ a/i - u4i o -o; 3 i 
1 

v a/ 1 - y 



21 



(36) 



0) 



This case is particularly meaningful, since it realizes one of the conditions (Q 
to saturate the bound (1231) for Si that is given by E\ = e(M 1 )Y 3 m atm /(_ft' 1 m^), with 
Y3 = — Imfco'ljJ [7]. Moreover the decay parameter is given by 



K\ — -f^min 1 1 — <^31 I + K & tvn I ^ 



L atm |^31 1 ? 



(37) 



where K atm = m atm /m* and K n 
then specialized as 



mi/m*. The expression (fl9|) for the projector gets 



mi |f4i| 2 |1 - lu^I +m 3 \U a3 \ 2 \ujl 1 \ + 2 y /m 1 m 3 Re[[/ al C/* 3 a/1 - 



a; 



31J 



la 



mi 1- + m 3 \<jj- 



(3? 



y 31l 



while, specializing the Eq. ([9]) for £i a and neglecting the term oc Xa , one obtains 



^atm 



i^i m 

^atm 

Kim,, 



U a 3\ 2 H jT - (|K*3| 2 — [ C^otl I 



m 



atm 



(39) 



mi m 3 



^atm ^atm 



mi + m 3 



m atm 



^31 



+ 



m 3 — m x 

^■atm 



Re 



Im 



^31 \/ 1 - U3I 
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31 



where m 3 /m atm = a/1 + m\jm 



''atm' 



If we first consider the case of fully hierarchical light neutrinos, mi = 0, then 



p _ £ la 
£ 1 



U, 



a3 I 



and 



2ei 



0. 



(40) 



Neglecting the effect of the running of neutrino parameters from high energy to low 
energy [27], one can assume that the U matrix can be identified with the PMNS matrix 
as measured by neutrino mixing experiments. We will adopt the parametrization [28] 



/ 



U 



Cl2 C13 



Sl2 C13 



" s i2 C23 — c i2 S23 Si 3 e tS C12 C23 — S12 S23 S13 e tS 



sise- 15 \ 

$23 c 13 



\ S12 s 2 3 - C12 c 23 S13 e % 



£2 



-C12 S 2 3 - S12 C 23 S13 e l C 23 C13 / 
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xdiag(e : 2 ; e 1 2 , 1) , 
(41) 



where = sin%, = cos 9ij and where we have used #13 = — 0.17, 6*12 = vr/6 and 
$23 = vr/4. One then finds 

P° < 0.03 , F% ~ J* ~ 1/2 . (42) 

This result implies a situation where the projector on the electron flavor is very small and 
the associated generated asymmetry as well, while the muon and tauon contributions are 
equal. Notice that even though there is a non-vanishing extra contribution to the muon 
and tauon CP asymmetries, they have equal absolute value but opposite sign. Therefore, 
since the projectors are equal as well, summing the kinetic equations fill I) over a, one 
simply recovers the one-flavor approximation with the wash-out reduced by a factor 2. 
This is a realization of the semi- democratic case that we were envisaging at the end of 
3.2 where K* ~ 7. In this situation flavor effects do not produce large modifications of 
the usual results, essentially a factor 2 reduction of the wash-out in the strong wash-out 
regime with a consequent equal relaxation of the lower bounds (see dashed lines in Fig. 
1). Moreover there is practically no difference between a calculation in the two or in the 
three-flavor regime. 

Let us now consider the effect of a non- vanishing but small lightest neutrino mass mi, 
for example m\ = 0.1 m atm . In this case the results can also depend on the Majorana and 
Dirac phases. We will show the results for io\ x purely imaginary, the second condition 
that maximizes the total CP asymmetry if mi <C m atm [7J. Notice that this is not in 
general the condition that maximizes r± a for non- vanishing mi, however we have checked 
that allowing for a real part of io\\ one obtains similar bounds within a factor 0(1) if 
mi < m atm . 

We have first considered the case of a real U. The results are only slightly sensitive to 
a variation of 9 is within the experimentally allowed range — 0.17. Therefore, we have 
set #i3 = 0, corresponding to U e3 = 0, in all the shown examples. In the left panel of Fig. 
4 we show the values of the projectors, the P° a 's, and of the ri a 's as a function of K\. 
The calculations are performed in the two-flavor regime since we obtain that successful 
leptogenesis is possible only for Mi > 10 9 GeV, where the two-flavor regime applies. Now 
a difference between the tauon and the sum of the muon and electron projectors and 
asymmetries arises. On the other hand for Ki ^> 100 this difference tends to vanish and 
the semi-democratic case is recovered again. 

In the central panel we have also plotted the quantities £i a (cf. (!26j) ) and their sum 
£1, where we recall that £1 gives the deviation of the total asymmetry from a calculation 
in the one-flavor approximation for hierarchical light neutrinos. One can see how now 
the contribution to the total asymmetry from the tauon flavor can be twice as large as 
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Figure 4: Dependence of different quantities on K\ for mi/m atm = 0.1 and real U. Left 
panel: projectors P® a and r lQ ; central panel: £i a and £1 as defined in Eq. fl26l) for 
thermal (thin) and vanishing (thick) initial abundance; right panel: lower bound on Mi 
for thermal (thin solid) and vanishing (thick solid) abundance compared with the one- 
flavor approximation result (dash-dotted line). 

from the electron plus muon flavor. For K\ ^> 100 the semi-democratic case is restored, 
the two contributions tend to be equal to the one-flavor approximation case and the total 
final asymmetry is about twice larger. Finally in the right panel we have plotted the 
lower bound on M\ and compared them with the results in the one-flavor approximation 
(dash-dotted line). At K± ~ 14 the relaxation is maximum, a factor ~ 3. For K\ ^> 
the relaxation is reduced to a factor 2, as in the semi-democratic case. 

Let us now study the effect of switching on phases in the U matrix, again for mi = 
0.1m atm . The most important effect arises from one of the two Majorana phases, $i. In 
Fig. 5 we have then again plotted, in three panels, the same quantities as in Fig. 4 for 
$i = — 7r. One can see how this further increases the difference between the e + /i and the 
r contributions and further relaxes the lower bound on Mi. The effect is small for the 
considered value mi/m atm = 0.1. However, considering a much larger value while keeping 
$i = — 7r, the effect becomes dramatically bigger. In Fig. 6 we have plotted the same 
quantities as in Fig. 4 and Fig. 5 for rn.i/m atm = 10. In the left panel one can see that 
now, for Ki ^> K min , |r le+/1 | ~ |r lr | ~ 2, much larger than in the previous case. This 
means that the dominant contribution to the flavored CP asymmetries comes now from 
the APic term. At the same time, very importantly, P° r P\ e +^ an d i n this way, as 
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Figure 6: Same quantities as in the previous two figures but with mi/m atm = 10 and o: 
non-vanishing Majorana phase: $i = — ir. 
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10" 1 10° 10 1 10 2 10 3 

Figure 7: Lower bound on M\. The solid lines are for vanishing phase, while the dashed 
lines are for $i = — ir. We have shown both the case of vanishing (thick lines) and of 
thermal (thin lines) initial abundance. 

one can see in the central panel, the dominant contribution to £i is given by £i T . This 
case thus finally realizes a one-flavor dominance case. The final effect is that the lower 
bound on Mi is about one and half order of magnitude relaxed compared to the one-flavor 
approximation. 

In the left panel of Fig. 7 we have summarized the dependence on mi/m atm of the 
lower bound on Mi, plotting both the case with zero Majorana phase and the case with 
$! = — 7i. One can notice how the the effect of the phase in relaxing the lower bound 
increases with mi/ro a t m . 

Finally we wanted to study the interesting case of a real orthogonal matrix Q, implying 
si = 0. If mi = there is practically no produced asymmetry, since the ei a oc E\ = 0. For 
a non-vanishing mi and a non-real U however the APi a ^ and consequently 6i a ^ 0. In 
Fig. 8 we have shown the results for mi/m atm = 0.1 and $! = ir/2. We have again plotted 
the same quantities as in Fig. 4, 5 and 6 in three different panels. One can see that an 
asymmetry can be still produced, as emphasized in [12] . However, successful leptogenesis 
is possible almost only in the weak wash-out regime. In the strong wash-out regime, for 
values Mi < 10 12 GeV, there is a very small allowed region only for K* ~ 14 < Ki < 30. 
This somehow means that, at least in the example we are considering and for small mi, 
flavor effects give a sub-dominant effect. 

We can thus draw two remarkable conclusions. The first is that the relaxation of the 
lower bounds compared to the one-flavor approximation is only 0(1) for mi m atm 
but grows up to one order of magnitude and more for mi > 10m atm ~ 0.5 eV, in the 
quasi-degenerate case. Therefore, when flavor effects are taken into account, one obtains 
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Figure 8: Same quantities as in Fig.'s 4,5 and 6 in the case of real fl for mi/m atm = 0.1 
and $i = 7r/2. The dot-dashed still refers to the unflavored case for Q = i? 13 and purely 
imaginary uj^. 

an opposite result compared to the one-flavor approximation where one has a suppression 
of the final asymmetry for growing absolute neutrino mass scale, leading to a stringent 
upper bound mi < 0.1 eV [291 111 EU E]- This upper bound seems now to hold only for 
Mi > 10 12 GeV [121 [13]. Notice however that a conclusive answer to this issue requires a 
full quantum kinetic calculation describing the evolution of the density matrices of leptons 
and anti-leptons [15]. The second is that, when flavor effects are included, leptogenesis 
is an interesting example of phenomenology, beyond /3/30v, where Majorana phases can 
play a relevant role. It must be said however that if the stringent cosmological upper 
bounds on the absolute neutrino mass scale will be confirmed, mi < 0.1 — 0.2 eV [30] . 
then flavor effects can only produce a 0(1) relaxation of the lower bound on Mi at large 
values of K\, compared to the results from the one-flavor approximation. This conclusion 
rigorously applies to the considered example but we arrived to the same conclusions also in 
other cases and it seems quite reasonable that this result holds on more general grounds. 
The reason is that for small mi <C m atm it does not seem possible to realize a one-flavor 
dominated scenario. Further studies are however needed on this point. 
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5 Final discussion 



Flavor effects in leptogenesis can produce significant deviations from the usual results 
obtained within the one-flavor approximation. However, these are quite restricted in the 
strong wash-out regime, imposing that the final asymmetry is unambiguously determined 
for different choices of the initial conditions. 

The traditional lowest bounds on M\ and T re h are essentially the same as in the one- 
flavor approximation. Therefore, flavor effects do not help to solve the conflict with the 
upper bound on the reheating temperature coming from the gravitino problem. For this 
purpose it is necessary to go beyond a hierarchical heavy neutrino spectrum [31], HJ [32], [21] . 

On the other hand, for a fixed value of the decay parameter K\, the lower bounds can 
be relaxed if it is possible to explain the observed asymmetry as dominantly produced in 
one flavor. This scenario seems not easily achievable since typically for low values of the 
projector, the final associated asymmetry is also lower. 

An interesting feature of flavor effects is that they depend in general on the low energy 
parameters [HI [12], in particular on the lightest neutrino mass m\ and on the phases in the 
neutrino mixing parameters. We have studied an interesting example assuming the condi- 
tions for maximal total CP asymmetry. We showed that for vanishing phases and lightest 
neutrino mass one essentially recovers the one-flavor approximation just with halved wash- 
out, the semi-democratic case, where the electron flavor gives a negligible contribution 
while the final asymmetry is equally shared between the muon and the tauon flavor con- 
tributions. However, switching on simultaneously a non-vanishing lightest neutrino mass 
and Majorana phases, lead to an enhancement of the CP asymmetry in the tauon flavor. 
In this way a tauon-flavor dominated picture is realized, where the wash-out is much 
lower compared to the case when flavor effects are neglected. It is then quite interesting 
that there is an emerging interplay between future low energy experiments and flavored 
leptogenesis predictions on the matter-antimatter asymmetry of the Universe. 
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Note added 

During the revision of the paper, the work [15] was completed, finding a more restrictive 
condition for the importance of flavor effects and for the validity of the Eq.'s ( fTTj) . The 
condition is particularly restrictive in the case of one-flavor dominance, relevant for the 
results presented in Figures 6 and 7 and for the possibility of circumventing the upper 
bound on the neutrino masses holding in the unflavored case. We should then warn the 
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reader that these results need to be confirmed by a more general quantum kinetic analysis. 



Appendix 

In this appendix we discuss the effects of AL = 1 scatterings [TFl [TEJ [19] on the deter- 
mination of the efficiency factors K la and more particularly on K±. It has been recently 
noticed that scatterings also contribute to the CP violating source term with the same 
CP asymmetries Si a associated to inverse decays (TQl EHJ, [14] . Therefore, here we also want 
to discuss the effect of this term extending the analytic procedure of [5] , where it was not 
included. Including scatterings the kinetic equations (fill in the iVxDS get modified as 

^ = -(D 1 + S 1 )(N Nl -N^ i ) (43) 

^ = e la (D 1 + St) (N Nl - N%) - Pi {W? + Wf) N Aa . (44) 

The notation is the same as in [5] , except that here we are keeping the index referring to 
the lightest RH neutrino. Defining 

JiW = l + ^, (45) 

the expression ( 1271) for n la in the case of initial thermal abundance (iVJy = 1) gets 
generalized as 

K{K la ) = fl-e ^ . 46 

ZB{Ki a )Ki a ]i(z B ) V / 

Comparing with the result obtained when scatterings are neglected in the CP violating 
term, there is now Ji{zb) instead of jf(zB), making the effect of scatterings even less 
important than it was already. This could be relevant in the case of a dynamically 
generated abundance in the weak wash-out regime, since scatterings are able to produce 
a larger abundance at the decay time. However again, when the effect of scatterings is 
included in the CP violating term as well, this enhancing effect is greatly reduced. The 
Eq. (!3~l!) for k_, when scatterings are included, becomes indeed 

Kl( Kl ,P? a ) . * e"^ ( e4^(*0 _ l) , (47) 



where 



SnK la t f 00 j / tt/id Si 



Ns(K la ) = —f^ + ^ dz> W£ g. (4f 
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The expression (|33|) for K f , remains formally unchanged, but now Ns(K\) has to be 
calculated replacing N(Ki) with N s (Ki) into the Eq. (|32|) . Even without specifying at 
all the scattering term Si, there is quite an interesting general feature resulting from the 
inclusion of scatterings in the CP violating term. One can indeed calculate the behavior 
of K la = k, + k_ for K\ — > discovering that [P 



cx Kf , (49) 



exactly as in the case when scatterings are neglected [5], except for the modified N Nl (z cq ). 
The reason is that now there is a proportionality between asymmetry generation and neu- 
trino production both in scatterings and in inverse decays. Conversely, when scatterings 
are neglected in the CP violating term, there is an overproduction of neutrinos due to 
scatterings not compensated by an equal (negative) production of asymmetry, such that 
scatterings enhance the final asymmetry that is oc K\ instead of K±. Notice however that 
in the strong wash-out regime there is practically no difference. 

In Fig. 9 we have plotted the final efficiency factor K la using a useful analytic expres- 
sion for ji (z) given by [5] 



z , / a\ Kf 
-In [! + -) + 1 



a V z J K\z 



1 + ^), (50) 



■ z 



where a = Kx/Kf ln(Mi/M^) and Kf ~ 0.1 K\. This approximation works well when 
scatterings involving gauge bosons and thermal effects [I8l [T9] are neglected. The different 
curves are obtained for Mi/Mjj = 10, 10 5 and 10 10 and compared with the case when 
scatterings are neglected (solid lines). Notice that the chosen values for Mi/Mg are the 
same as in [5] . A comparison shows that the results exhibit now a greatly reduced model 
dependence, here illustrated by the different values of Mi/M^, also in the weak wash- 
out regime. Including thermal effects and scatterings involving gauge bosons would yield 
similar results. 

In conclusion, the inclusion of scatterings in the CP violating terms greatly reduces the 
theoretical uncertainties in the weak wash-out regime, even though these are still much 
larger than in the strong wash-out regime. In particular the asymptotic behavior oc K 2 at 
small K is a robust feature, as explained by the Eq. (j4"9l) . In the right panel we have shown 
the corresponding value of K±. Notice that thermal effects are approximately reproduced 
by a thermal Higgs mass Mh — 0.4 T [THl HE] and in the strong wash-out regime this 
is approximately equivalent to the case M 1 /M H ~ 10 where scatterings produce just a 
correction in the calculation of K+. Even for larger M\jMn the value of K+ calculated 
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Figure 9: Efficiency factor (left panel) and K* (right panel) when scatterings are taken 
into account for 3 values of Mi/Mr: 10 (dash), 10 5 (dot-dash), 10 10 (dot-dot-dash), to be 
compared with the case when they are neglected (solid lines). 



neglecting scatterings is within the precision needed for the results of the paper, thus 
justifying the approximation to neglect scatterings. 
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